Electron optics with magnetic vector potential barriers in graphene 
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An analysis of electron transport in graphene is presented in the presence of various arrangement 
of delta-function like magnetic barriers. The motion through one such barrier gives an unusual 
non specular refraction leading to asymmetric transmission. The symmetry is restored by putting 
two such barriers in opposite direction side by side. Periodic arrangements of such barriers can be 
used as Bragg reflectors whose reflectivity has been calculated using a transfer matrix formalism. 
Such Bragg reflectors can be used to make resonant cavities. We also analyze the associated band 
structure for the case of infinite periodic structures. 

PACS numbers: 81.05.Tp,73.23.-b,73.63.-b,78.20.Ci 



i 

C 

o 
o 



^1- 

> 

in 

On 
(N 

o 
oo 
o 



X 



In a two-dimensional electron gas, there is a well es- 
tablished similarity between ballistic electron transport 
through electrostatic potential barriers and light propa- 
gation in geometrical optics [l| . This has been extended 
to the massless Dirac fermions in Graphene 0| where it 
was recently established that electron transport in the 
presence of an electrostatic potential barrier is analogous 
to negative refraction through metamaterials0, 0|- The 
relativistic behavior of Graphene electrons also leads to 
Klein tunneling [H, [fjj], where a relativistic particle can 
tunnel through a high barrier by the process of pair pro- 
duction and thus it is not possible to confine them using 
such potential barriers. Such confinement is however pos- 
sible by using magnetic barriers [5~oj | . Can we understand 
this difference in behavior of massless Dirac fermions by 
comparing to propagation of light? The difficulty in us- 
ing an optical analogy is that, unlike the electrostatic 
potential, the magnetic vector potential couples with the 
momentum of the electron. 

In this work, we show that wave- vector dependent tun- 
neling of massless Dirac fermions through magnetic barri- 
ers 0, EH can be understood in terms of well-known 
ideas in geometrical optics. However, the correspond- 
ing SnclPs laws are very different from those of ordinary 
geometrical optics. We then carry out this analysis to 
propose devices such as a Bragg reflector using a transfer 
matrix approach and qualitatively depict how a resonant 
cavity can be constructed with such a reflector. Addi- 
tionally, we comment on the band structure of electron 
transport when such magnetic barriers are placed peri- 
odically. 

The proposed structure consists of a graphene sheet 
placed in close proximity to long magnetic stripes that 
produce highly localized magnetic fields as depicted in 
Figffl Such field profiles can be created using demag- 
netizing fields produced at the edges of narrow stripes 
made with hard ferromagnetic materials of perpendic- 
ular or in-plane anisotropy. It is possible to make such 
stripes at various length scales. Materials such as CoCrPt 
used in magnetic recording produce field strengths of 1 
Tesla close to the surface with bit lengths ranging from 



50-100nm. Patterned stripes down to lOnm can be real- 
ized using nanolithography It is possible to achieve 
field profiles at even smaller dimensions using domain 
walls with widths in the range of 10-50nm and magnetic 
nanostructures down to 5nm [3, 14 1 and 0.15nm [HI 



having highly localized field variations. The particular 




FIG. 1: Monolayer graphene with ferromagnetic stripes hav- 
ing magnetizations perpendicular (a) and parallel (b) to plane. 
The magnetic field B (red) produces a magnetic vector poten- 
tial A (blue). Single MVP barriers are formed in (a) and 
(b). Also shown are a double MVP barrier (c) and a periodic 
lattice (d). 

form of magnetic barrier that we shall use in our calcu- 
lation can be realized using two narrow ferromagnetic 
stripes of perpendicular anisotropy with appropriately 
narrow dimensions and magnetized in opposite directions 
(FigQJa)). The same profile can also be achieved by one 
ferromagnetic stripe whose magnetization is parallel to 
the graphene sheet at a height zq above it (FigJTJb)). 
Such barriers have been used in the literature 0, @] and 
the magnetic field of such structures is 

B = B(x, zq)z = Bq[K(x + d, zq) — K(x — d, zq)]z 

Here, K(x,Zq) = — and Bq is a constant depen- 

dent on the aspect ratio of the stripe. For a given value 
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of zq we plot in FigQJa) the profile of such a magnetic 
field and the corresponding vector potential. As the plot 
shows, such an inhomogenous magnetic field can be well 
approximated as a delta function-like magnetic barrier. 
This approximation is valid for all the length scales dis- 
cussed in the preceding paragraph as long as the typical 
magnetic length l B = \J is much larger than the 
width of such magnetic barriers. Accordingly, we use 
such delta function-like barriers for the rest of the paper. 
This choice is guided by the fact that it is amenable to 
simpler mathematical treatment, thus making the anal- 
ogy more transparent. We shall, however, point out later 
that the discussed analogy with geometrical optics is 
more general and is applicable for other types of mag- 
netic barriers with finite width I IL Oil . 



Here, A = 1 for |sc| < d and = for |sc| > d. The above 
two coupled equations can be decoupled easily and result 
in a Schrodinger like equation of the form 

[-^ + (y B - A) 2 ]0!, 2 = e 2 0i, 2 . 

In the region — d < x < d electrons see a barrier of height 
[fcj,+sgn(e) j^] 2 ■ The corresponding wavefunctions in any 
region of space can be written in terms of a linear super- 
position of forward and backward moving plane waves 
such that 

! e ik x x _|_ re ~ik I x x ^ 
ae iq*x + be -iq*x \ x \ < d (3) 
te ik„x x > ^ 
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FIG. 2: (Color online) (a) Asymmetric refraction through a 
single barrier, (b) Refraction through a double barrier where 
the symmetry in the trasmission is restored. Barrier regions 
which are denser or rarer in terms of refraction of the wave 
vector are shaded with the same color. 

We use the following magnetic field and vector poten- 
tial in Landau gauge for a magnetic potential barrier Q : 



B = 



,[S(x + d)-S(x-d)]z;A v (x) = B£ B Q(d 2 -x 2 )y 



Since Q(x) is the Heaviside step function, we call this a 
magnetic vector potential (MVP) barrier. 

For massless Dirac fermions in graphene in MVP bar- 
riers, we consider the limit where the electrons at the 
K and K' points are decoupled from each other [lij ]. 
Near each such point, the wavefunction is given by a 
two-component spinor and satisfies the equation 



v f (tt x ± in v )tp 



2,1 



1.2 



(1) 



Here, vp is the fermi velocity (« c/300) and 7r = p+ ~A 
Using as the unit of energy such that e = ^jf, £p 
as the unit of the length such that x = -f^ , sgn(e) = — 1 
and i/j = 4>{x)e lky v in the Landau gauge, we get 



d 

i[— ± (k y t B ~ A)] 02,1 = £01,2 



(2) 



s [ e i(Kx+<t>) _ re -i(fex+«] x <-d 
. s '[ ae i{i^+6) - be-i^x+e-jj \ x \ <( i (4) 
ste i(Kx+<t>) x> d 

Solutions of the above equations arc very different from 
those in the presence of a uniform magnetic field since 
here the magnetic field is highly non-uniform and has 
singular delta function like structures. Also, s, s' are 
given by sgn(e) and are both +1 for electrons when only 
magnetic fields are present and no electrostatic poten- 
tials are applied. A similar treatment can be done in 
the presence of an additional electrostatic potential, in 
which case both s,s' are not necessarily +1 [l7j . The 
wave vector components are [k x ,k y ] = fcp[cos</>, sirup] 
outside the magnetic barrier and is the incident an- 
gle for an electron wave. The fermi energy of the in- 
cident electrons is Ep = hvpkp. In the dimensionlcss 
form, this is ep = kp£ B and will control the trans- 
port by changing the refractive index of the barrier re- 
gion. Since magnetic field does not do any work, en- 



ergy conservation gives fc 2 , 
for \x\ < d kp sin ( 



k 2 



kp for 



-) 2 = kp for \x\ < d. Since 



|ar| > d and 

COS' 



-l(3»_) 



sin \6\ — sin |0| — sgn(0) 



gives 



1 



7r^'-?<^<J (5) 

kpis 2 2 



The situation has been depicted in Fig. [Ha). The wave 
incident with positive wave vector will bend towards 
the normal whereas the waves incident with negative in- 
cidence angle, corresponding wave vector will bend away 
from the surface normal inside the barrier region. There- 
fore, the Snell's law of electron waves in such magnetic 
barriers is not specular as it is for light wave on smooth 
surface or for the incidence of the electrons on an elec- 
trostatic potential barrier [3, HJ . This unusual refraction 
can also be thought as a consequence of breaking of time 
reversal symmetry in presence of such magnetic barriers. 
When the magnetic field will be reversed, the denser and 
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rarer medium will interchange their side without chang- 
ing the asymmtric transmission. 

According to Eq. [5l when |sin|0|| > 1, becomes 
imaginary and the wave in the second medium becomes 
evanescent. In the language of optics this corresponds to 
total internal reflection (TIR). According to Fig. [Ha), 
this will happen when sin|0| > 1 for — $ < <j> < and 
when sin 1 1 < — 1 for < < -|. In the latter case, 
this requires the wave vector to be negatively refracted 
3 at sufficiently high magnetic field before TIR occurs. 
It also follows that for a given strength B the magnitude 
of critical angle of incidence |0| = C for TIR is higher 



for < < ? as compared to the one for 



< < 0. 



Because of TIR the transmission on both sides of Fig [3] 
drops to beyond a certain value of and this value is 
lower for negative angles of incidence . 

The wavefunctions given in Eqs. [3] and [4] are similar 
to those of massless Dirac electrons scattered by an elec- 
trostatic step potential considered in [H| . This is because 
the MVP barrier creates a momentum dependent step 
potential of [k y + sgn(e)^] 2 . Continuity of the wave- 
function at the boundaries of the MVP barrier can be 
used to calculate the transmission coefficient as 

2ss'e~ lkxD coscbcosO 
t = 

ss >[ e -iq*D cos (0 + 0) + e tq*D cog (0 _ 0)] _ 2i smq x D 

(6) 

Here, D = 2d. Thus t, transmittance T = t*t and re- 
flectance R = 1 — T have same expressions for electro- 
static potentials in Ref. There are, however, key dif- 
ferences. For electrostatic barrier as — > — 0, 9 — > —6. 
Here because of EqJS] that is not the case. Thus, the 
same Eq|6] gives symmetric transmission for electrostatic 
potential in [f| and asymmetric transmission here. 
For high electrostatic barriers such that V ^> Ep, the 

wave vector is given by q x = ^ ^f^Jf - k 2 , which is 

real. The corresponding transmittance for electrostatic 
potentials is 



T = 



cos" 



1 — cos( q x D) sin 



(7) 



and is 1 at — 0. This exhibits Klein tunneling for 
massless Dirac fermions For the magnetic barrier, 

, \ tx y/B. However, unlike the electrostatic field case, 
F B . . . 

the magnetic field changes the wave vector in both di- 
rection and not the energy. At high magnetic field, 
1x = hp - (k y — j^) 2 = — k 2 < 0. As discussed, this 
leads to TIR and not Klein tunneling. In Fig. [3J the 
magnitude of critical angle beyond which TIR occurs is 
lower for a higher magnetic field. Then, a stronger MVP 
barrier leads to higher reflections as opposed to complete 
transmission at normal incidence by high electrostatic 
potential barrier. A similar situation is encountered in 
other forms of magnetic barriers [TFJ, [TJI • 



Complete transmission only occurs for q x D = m: in 
Eqj6] This corresponds to resonant tunneling for Dirac 
electrons and happens in the same way as for non- 
relativistic electrons, appearing as a number of peaks in 
the plots in Fig[3J The number of such tunneling peaks 
increases with barrier width for both MVP barriers and 
electrostatic barriers. 




FIG. 3: Polar plot of T vs. <f> for single MVP barriers 
of widths (a) lOOnm, (b) 500nm. Blue (thick) =0. IT with 
lB=81nm, Black (thin)=3T with /s=14.8nm. 

To get symmetric transmission out of such a barrier it 
is therefore required to place two such MVP barriers side 
by side but oppositely oriented as depicted in FigQTc). 
The magnetic field creating such a barrier is 



B = B z (x)z 



,[6{x + d)+ S(x -d)- 2B6(x)]z 



We again consider energy conservation in medium 1 
(—d < x < 0) and medium 2 (0 < x < d), which gives 



Ql,2 



fc F ;sin|0i, 2 | 



sin |0| Tsgn(0) 



1 



kpii 



(8) 



The incident angle is 



< <t> < 5 and the angle of 



refraction is 6\ and 62 in media 1 and 2 respectively. The 
absolute value of relative refractive index of region 1 with 
respect to region 2 on the left side of the surface normal 
is just the inverse of that on the right side of the surface 
normal and can be combined in the following expression 



sin 



1"2 = 



sin |0| - sgn(0) 



k F l t 



sin 1 6» 2 1 sin |0| +sgn(0)^- 



(9) 



Thus, for such double MVP barriers, the wave vector 
bending towards (away from) the surface normal in the 
first half of the barrier bends away from (towards) the 
surface normal in the second half of the barrier as shown 
in Fig[2jb) . This will achieve symmetric transmission 
through such a barrier as demonstrated in Fig[5] Criti- 
cal angles of incidence beyond which TIR will occur for 
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positive and negative <f> will also be interchanged as we 
go from the first to the second barrier region. However, 
at higher B fields, TIR will occur at both regions of the 
barrier. Consequently, the total reflectivity of the barrier 
increases as can be seen by comparing Figs [3] and HI 



Here / is the unit matrix and 



Mi 



(12) 



The solution of these continuity equations can then be 
written as 




FIG. 4: Polar plot of T vs. <f> for DMVP barriers of widths (a) 
lOOnm, (b) 500nm. Blue (thick)=0.1T with Ib =81nm, Black 
(thin)=3T with Zs=14.8nm. 

Practical devices such as Bragg reflectors can be made 
by exploiting the high reflectivity of double MVP barriers 
to manipulate electrons, as will be discussed later. Such 
structures, if realized, could also be very useful to confine 
carriers in desired areas in graphenc and away from edges, 
where edge states could adversely affect transport. 

To calculate transmittance for double MVP (DMVP) 
barriers, we write the wavefunction in the same way as 
in Eqsl5]and|U 



_|_ re ik*x x < _^ 

ae lqix + be~ lqix x G [-d, 0] 

ce iq2X + de~ lq2X x G [0, d] 

te %k,cX x > d 



(10) 



s ^ e i(k x x+if>) _ re -i(fe x a;+0)j x < - 

s 1 [ae i ^ lX+e ^ - be~ i{ - qxX+e ^] x G 



S2 ce 



i(q 2 x+9 2 ) 



-d 

-d,0] 

ste i(k I x+<t>) x> d 



(11) 



The corresponding sign factors associated with <f> 2 in 
these regions are s\ and s 2 and are both 1. The transmit- 
tance and reflectance can now be easily computed by im- 
posing the continuity conditions on the above functions 
at the location of the barriers, namely at x = —d, 0, d. We 
can express the transmittance and reflectance in a com- 



pact form by introducing A = e 
to define the following matrices: 



and B\ ; 



-iqi t2 x 



M A 



A 
A* 



1 



M s 



1 

s,si,s 2 



, M Sl 2 ^ 

= I, M Bl>2 = 



1 

-i@l, 2 ,4> 

B h2 
BU 



[1 r] T = Mj^MtjT 1 [M$ 1 Mb-l Mg^ 

M^Mb^M^M^Ma* [t 0] T 

To understand the above formulae, we introduce the 
transfer matrix through a DMVP barrier: 

Tdmvp = M Sl M Bl M- 1 Me 2 M B2 M- 1 

We can interpret M7 1 Mg 1 as the phase shift at the first 
boundary, Mg^M^ as the phase shift at the last bound- 
ary and Mg i 1 M$ 2 as the phase shift at the barrier at 
x = 0. Thus, we can rewrite the above equation as 



[1 rf = MZM^TbmvpM+Ma. [t 0] 



(13) 



The resulting two equations can be solved to yield the 
transmittance through DMVP barriers, which has been 
plotted in FiglH 

We shall now discuss similarities between transport of 
masslcss Dirac fcrmions through MVP barriers with elec- 
tromagnetic propagation in periodic stratified media (lij 
as well as Dirac fermions in periodic electrostatic poten- 
tials [23]. In real structures, there will only be a finite 
number of barriers and the lattice translational symme- 
try will break down at the boundary. To simplify the 
analysis, we assume the DMVP barrier structure can be 
repeated infinitely We consider each unit cell of size 
D = 2d. The n-th cell is given by (n — l)D < x < nD. 
In the a-th part of the given unit cell, the wavefunction 
is 

4>x = a*e %q «* {x - nD) +6° e - i9 « (a, - n£,) 
02 = < [a°e*fo~ (*-*»)+*«] ~ bZe- i[q ^ ix - nD)+9a] 



n 



r h 2 — <t a n - 
inD rcvca J s tfl C 



Here, a = 1, 2, a* = a„, b\ = b n , a\ 
qi , = <? 2 • The exponential factor e 
lattice translational symmetry, which is not present for 
isolated single and double barrier structures. Imposing 
continuity at both interfaces of the n-th unit cell gives 



Ma^Mg 



Cn-1 

d n -i 



M 3ln M 9l M Bl 



M Sln M ei M E 



M S2 M 02 M B2 



a r , 
bp 

d„ 



(14) 



Imposing Bloch condition gives the band structure from 
the following eigenequations 





= e} KD 


Cn-1 




~K n 


K 12 ~ 




Cn-1 


d n _ 




d n -l 




K 2 i 


K 2 2_ 




d n -\ 
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where K is the Bloch momentum and the matrix ele- 
ments Kij can be calculated from Eg [Ml as 



K n K12 
K21 K 2 2 



= (Mg 2 Mb 2 ) ~ 1 Mg 1 {Mq 1 Mb % ) ~ 1 Mg 2 



This is the same as T D ^vp permuted 
det Ki 



Unitarity gives 
1, yielding the eigenvalue equation 

KfaB) = i cos" 1 [^TriKij)] 

The condition \^Tr(Kij)\ < 1 corresponds to propagat- 
ing Bloch waves whereas |iTr(ify)| > 1 leads to evanes- 
cent Bloch waves that correspond to forbidden zones in 
the band structure. Such band structures have previ- 
ously been studied for many problems including con- 
densed matter systems, optics [19j and relativistic quarks 
[20I ]. While the detailed band structure is provided else- 
where [17|, we plot in FigO ||Tr(ify)| as a function of 
the incident angle </> for different B values. A forbidden 
region appears at <f> — at higher B due to larger differ- 
ence between the refractive indices of adjacent regions. 



B=0.1T 




-7r/2 



-tt/4 tt/4 

B=3T 



jr/2 



waves similar to the ones given earlier in Eqs. [TOl and [TT1 
for one DMVP barrier. To solve this set of equations, we 
proceed as earlier using continuity of the wave function 
at the magnetic barriers at x = pd, —n < p < n. Just as 
Eq[T2] describes the solution for one DMVP barrier, the 
solutions for n DMVP barriers can be written in matrix 
form as 



[1 



{M A - l TM- l TZ MYP M^M A *Y[t 0] T (16) 



Here, we have used Eq ll2l and Tdmvp is just the transfer 
matrix through a DMVP barrier. We can again inter- 
pret MT Mg x as the phase shift at the first boundary 

and MqM,), as the phase shift at the last boundary. 
These two boundaries are special since these are the in- 
terfaces of the magnetic medium with the non-magnetic 
region. What appears in the middle is the transmission 
through DMVP barriers repeated n times. A represen- 
tative plot is given in Fig[5J As can be seen, a practical 
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FIG. 6: Reflectance versus (j> through a Bragg reflector with 
different periods of DMVP barriers of d—100nm and B~0.1T. 



FIG. 5: cos(A'D) vs. <f> for an infinite periodic lattice 
with d—100nm showing evanescent (dashed) and propagating 
(solid) Bloch waves. For comparison with earlier plots, cj> is 
still used. 

We shall now modify the above result of infinite peri- 
odic barriers to analyze a finite chain of DMVP barriers 
that makes a Bragg reflector. A Bragg reflector can be 
formed by superposing side by side n such DMVP bar- 
riers. We shall briefly describe how the transmittancc 
and reflectance of such a reflector can be calculated us- 
ing a transfer matrix formalism. The magnetic field for 
a Bragg reflector placed symmetrically around the origin 
can be written as 



B = B z (x)z 



?b[8(x + nd) +6(x- nd) 

n-l 



+ {-l) p+n 2BS{x - pd))]z (15) 

p— 1 — n 

The series of wave function solutions in the various re- 
gions are linear combinations of right and left moving 



Bragg reflector with a high enough reflectance can be re- 
alized with just a few periods of DMVP barriers. A Bragg 
reflector with large n is broadly similar to an infinite pe- 
riodic lattice. Particularly at low B (0.1T) around = 0, 
there is high transmission and R = for both structures. 
Similarly, at high B (3T), there is strong supression of 
transmission near d> = in both cases. 




FIG. 7: Variation of current (T(B)) through a Bragg reflector 
with period n. B=Q.1T (A), IT (*), 3T (O) and d=100nm. 



6 



The higher reflectance of a Bragg reflector will strongly 
suppress transport. We define the average transmission 
as transmittance T(<J)) multiplied by the cc-component of 
velocity integrated over all angles of incidence for a given 
B and d such that 

(T(B)) =2v F [* d<t> cos <j>T{<j>) (17) 
Jo 

FigE] plots (T(B)) for various B and d and shows trans- 
mission is strongly suppressed with increased magnetic 
field and this happens within a few periods of the Bragg 
reflector. The above formula, when generalized for a 
range of energy levels, leads to the conductance of the 
structure, which could be measured experimentally. 

We have thus shown that reflectance can be controlled 
by suitably modifying the strength and locations of the 
magnetic barriers and thereby changing the refractive in- 
dex of the intervening medium in a novel manner. This 
principle could be the basis of more elaborate structures 
depicted in FiglU In a magnetic waveguide (Fig[8ja)), 
reflection must be high at desired propagation angles, 
which could be manipulated by changing the magnetic 
field. For the resonant cavity shown in FigJSfb), high 
reflection is needed near normal incidence. Geometries 
such as three-mirror or four-mirror cavities could be used 
for high reflection at other angles. 




FIG. 8: Bragg reflectors with MVP barriers used as a mag- 
netic waveguide (a), and as resonant cavities (b-d). 

To conclude, we have shown that transport through 
MVP barriers can be understood in terms of propagation 
of light through periodic stratified media. This analogy 
can partially be attributed to the fact that Eqj2] describ- 
ing Dirac fermions and the Maxwell equations [2lJ are 
both linear wave equations, although there are some key 
differences that have been pointed out in this paper. 

This formalism describes transport in the ballistic 
regime, which corresponds to the case of pristine, low- 
doped Graphcne. The same treatment can be also ex- 
tended to non-relativistic electrons. This picture will 
be modified when the effects of disorder and electron- 
electron interactions arc included. 



Using these concepts, practical devices such as Bragg 
reflectors for manipulating Dirac electrons in graphcne 
can be made. Such barriers suppress Klein tunneling, 
thereby achieving confinement in graphcne which can be 
seen through strong supression of transmission of elec- 
trons. 
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